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ABSTRACT 

In this report, the heat transfer and isothermal flow friction 
characteristics for seven compact heat exchanger surfaces are presented 
and compared. Five of the surfaces were of a similar triangular fin 
configuration, differing in material and method of construction, and 
two of the surfaces, having been tested previously, were of a modified 
rectangular passage configuration. The modified rectangular passage 
surfaces showed better overall performance than the triangular fin sur- 
faces, and matrices constructed from perforated material showed better 
overall performance than similar surfaces constructed from non-perforat- 
ed material. For similar matrix configurations, small variations in 
manufacturing procedure had little effect upon heat transfer and flow 
friction performance. 

The heat transfer data was obtained by means of the single-blow 
transient testing technique. A cylic method of transient testing was 
also investigated, and found to be both reliable and useful in that the 
additional method of testing extends the Re3molds number range for the 
experimental apparatus. 




TABLE OF CONTENTS 



&%Tas„r°°^ 



SECTION 




Page 


1. 


Introduction 


13 


2. 


Summary of Theory 


14 


3. 


Experimental Techniques 


23 


4. 


Descrlf>tlon of Matrices 


27 


5. 


Presentation of Results 


29 


6. 


Experimental Uncertainties 


31 


7. 


Discussion of Results 


33 


8. 


Recommendations for Further Work 


38 


9. 


Bibliography 


39 


APPENDIX 


A Data Reduction Relationships 


88 


APPENDIX 


B Description of Equipment 


95 



APPENDIX C Conduction Parameter 

APPENDIX D Calculation of Experimental Uncertainties 
APPENDIX E Computer Program for Data Reduction 



98 

100 

101 



3 



LIST OF ILLUSTRATIONS 



FIGURE 




PAGE 


1. 


Schumann's curves 


42 


2. 


Lockers curves 


43 


3. 


Locke's curves 


44 


4. 


Ntu versus Maximum Slope 


45 


5. 


Ntu versus Maximum Slope 


46 


6. 


Ntu error versus Ntu 


47 


7. 


Ntu versus Amplitude Ratio 


48 


8. 


Schematic of Test Apparatus 


49 


9. 


Photograph of Test Apparatus 


50 


10. 


Nichrome Wire Heater 


51 


11. 


Inlet, Heater, and Test Section 


52 


12. 


Test Section with Matrix Holder Removed 


53 


13. 


Matrix Holder 


54 


14. 


Apparatus for Cyclic Testing 


55 


15. 


Effects of Longitudinal Conduction Parameter 


56 


16. 


#1 Geometric and Physical properties 


57 


17. 


Solar #2 Geometric and Physical properties 


58 


18. 


S6'lar #3 Geometric and Physical properties 


59 


19. 


St>hi*r #4 Geometric and Physical properties 


60 


20. 


Solar # 5 Geometric and Physical properties 


61 


21. 


160/40 TV Modified Rectangular Passage 


62 


22. 


Nickel Modified Rectangular Passage 


63 


23. 


j and f versus N__ Solar ,'>1 
RE 


64 


24. 


j and f versus Solnr #2 

RE 


65 



4 



FIGURE 



PAGE 



25. 

26. 

27. 

28. 



j and f versus N_„ Solar 
RE 

j and f versus N,,„ Solar #4 
j and f versus N_„ Solar #5 

KJi 

j and f versus N 160/40 TV Modified Rectangular Passage 



66 

67 

68 
69 



29. j and f versus N Solid Nickel Modified Rectangular Passage 70 

30. Comparison of Heat Transfer Data for all Surfaces 71 

31. Comparison of Flow Friction Data for all Surfaces 72 

32. Comparison of Flow Area Goodness Factor for all Surfaces 73 

33. Comparison of Heat Transfer and Friction Power for all 74 

Surfaces 

34. Heat Transfer Data by Means of Cyclic Method Solar #3 75 

35. Heat Transfer Data by Means of Single-Blow and Steady 76 

State Steam to Air Test for Solar #2 

36. Comparison of Experimentally Determined Heat Transfer 77 

and Friction Data with Theory for Solar #3 

37. Chart Trace from Cyclic Test 78 

38. Conduction Parameter Analog 79 



5 



LIST OF TABLES 



TABLE PAGE 

I, Sumary of Heat Transfer and Flow Friction Results 80 

SOLAR #1 

II. Summary of Heat Transfer and Flow Friction Results 81 

SOLAR #2 

III. Summary of Heat Transfer and Flow Friction Results 82 

SOLAR #3 

IV. Summary of Heat Transfer and Flow Friction Results 83 

SOLAR #4 

V. Summary of Heat Transfer and Flow Friction Results 84 

SOLAR #5 

VI. Summary of Heat Transfer and Flow Friction Results 85 

160/40 TV Modified Rectangular Passage 

VII. Summary of Heat Transfer and Flow Friction Results 86 

Solid Nickel Modified Rectangular Passage 

VIII. Ntu as a Function of Maximum Slope and Longitudinal 87 

Conduction Parameter 



6 



NOMENCLATURE 



English Letter Symbols 




A 


Matrix total heat transfer area 


sq ft 




Base area (plane surface area) 


sq ft 


A 

c 


Matrix minimum free flow area 


sq ft 


Ai: 

fr 


Matrix total frontal area 


sq ft 


A 

o 


Open area of matrix material due to 
perforations 


sq ft 


A 

s 


Matrix solid cross sectional area 
available for thermal conduction 


sq ft 




Matrix solid cross sectional conduction 
area considering effect of matrix material 
perforations 


sq ft 


c 

p 


Specific heat (gas) at constant 
pressure 


Btu/(lbm deg F) 


c 

s 


Matrix material specific heat 


Btu/(lbm deg F) 




Flow stream capacity rate (roCp) 


Btu/(hr deg F) 


c 

s 


Matrix heat capacity 


Btu/(deg F) 




Hydraulic diameter of internal 
passage 


ft 


E 


Friction power expended per unit of 
surface heat transfer area 


hp/sq ft 


F 


Frequency of sinusoidal temperature 
oscillation 


cycles/sec 


G 


Matrix flow stream mass velocity (m/Ac) 


Ibm/ (hr sq ft) 


8c 


Proportionality factor in Newton’s 
second law 


32.2(lbm ft)/(lbf 



h Thermal convection surface heat transfer Btu/(hr sq ft deg F) 
coefficient; heat transfer power per 
unit area per degree temperature difference 

K Loss coefficient for sudden flow contrac- dimensionless 
^ tion at matrix entrance 
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K 

e 


Loss coefficient for sudden flow 
expansion at matrix exit 


dimensionless 




k 


Unit thermal conductivity 


Btu/ (hr sq ft 


deg F/ft) 


k 

s 


Matrix material thermal conductivity 


Btu/(hr sq ft 


deg F/ft) 


L 


Total matrix flow length 


ft 




m 


Mass flow rate 


Ibm/hr 




P 


Pressure 


Ibf/sq ft 




P 


Matrix porosity 


dimensionless 




q 


Heat transfer rate 


Btu/hr 




R 


Universal gas constant (53.35 for air) 


(ft lbf)/(lbm 


deg R) 


r 


Electrical resistance 


ohm 






Hydraulic radius (A^L/A^) 


ft 




t 


Temperature 


deg F 




u 


Flow velocity 


ft/sec 




V 

m 


Matrix volume 


cu ft 




V 

s 


Matrix material volume (corrected 
for effects of perforations) 


cu ft 




W 

s 


Matrix mass 


Ibm 




X 


Distance from matrix inlet in 
flow direction 


ft 
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Greek Letter Symbols 



A 

nr 

6 

/ 

V 

(® 

cr 



Area density; compactness (ratio of matrix total 
heat transfer area to core volume-A/V ) 

Difference or change (time, distance, temperature) 

Fin efficiency 

Time 

Fluid viscosity (dynamic) 

Open area ratio due to matrix material 
perforations 

Density 

Ratio of free flow to frontal area (A /A- ) 

c fr 

Denotes "function of" 

Angular frequency of temperature oscillation 



sq ft/cu ft 



Sec, hr. 
Ibm/hr ft 

Ibm/cu ft 
dimensionless 

rad/sec 
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Subscripts 



atm Local atmosphere 

ave Average 

£ Fluid (gas, air) 

i Initial, inlet, individual 

k EqullllKnt 

m Mean, matrix 

0 At orifice 

s Solid (matrix material) , static 

STD Standard (temperature and pressure) 

w Wall (solid surface) 

X Local conditions 

1 At inlet (upstream of heating elements) 

2 At matrix inlet (downstream of heating) 

3 At matrix exit 
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Dimensionless Groupings 



f 



j 



M 



N, 



Nu 



N 



Pr 






Ntu 



R. 



Fanning friction factor (Ratio of wall shear stress to the 
fluid d)niamic head)* This factor is plotted versus Nr to 
illustrate the flow friction characteristics of a matrix 
surface* 

'ij 

Colburn Factor ( Np^ ^ )> a generalized heat transfer 
grouping. This factor is plotted versus Nr to illustrate heat 
transfer characteristics of a matrix surface 



rr\ c, 



Flow parameter ( ' 

Nusselt number ( ) , a heat transfer modulus 

Prandtl number a fluid properties modulus 

Reynolds number ( ^ flow modulus 

Stanton number ( )» a heat transfer modulus 

hA/ 

/me 



Number of heat transfer units ( 
parameter 






) , a heat transfer 



/ 

Ratio of matrix outlet temperature amplitude to matrix inlet 
temperature amplitude during cyclic temperature variation 






X 

Xk 

n: 



Longitudinal heat conduction parameter for solid matrix 
material ( >• 

Longitudinal heat conduction parameter for perforated matrix 
material ( >• 

Distance parameter (x/L) 

Time parameter ( ) used in Locke's analysis 

Time parameter ( ^ ) used in Cyclic technique 

Frequency parameter ( ) used in cyclic technique 
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1 . 



INTRODUCTION 



A compact heat exchanger may be defined as a heat exchanger having 
a large heat transfer surface area per cubic foot of core volume* 
Arbitrarily, heat exchangers with more than 200 square feet of heat 
transfer area per cubic foot of core volume are considered compact* 
Applications of particular interest for compact heat exchangers, are gas 
turbine power plants in aircraft, automobiles, and small seagoing craft. 

In all these applications the designer is limited by weight, flow frontal 
area, core volume, and pressure drop across the heat exchanger|^ An effec- 
tive basis for presentation and comparison of heat transfer and flow 
friction data for various heat transfer surfaces and core configurations 
is by plotting the dimensionless groupings- Colburn's heat transfer fac- 
tor, "j", and the Fanning friction factor, "f“, versus Reynolds number for 
heat exchanger surfaces. 

The data presented in this report was obtained by means of the 
single-blow transient testing technique described in the two sections 
following, and in reference [24]* The objectives of this thesis were 
to obtain and compare flow friction and heat transfer data for compact 
heat exchanger matrices of different configuration, to correlate this 
data obtained with data obtained for identical surfaces by means of 
a steady state steam- to-air testing technique, and to develop a testing 
technique that will extend the range of transient testing to higher 
Reynolds numbers. 
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2 . 



SUMMARY OF THEORY 



A. BACKGROUND 

Analytical development of the theory for the single-blow or 
transient heating/cooling problem was performed in the 1920 *s by Anzelius 
(1926), Nusselt (1927), Hausen (1927-1929), and Schumann (1929). 

There was some duplication of effort in that Nusselt did not know 
of Anzelius* publication, and Schumann was unaware of the work of 
Anzelius and Nusselt. Anzelius and Schumann used similar procedures. 
Anzelius represented the gas and solid temperatures as unsolved, re- 
latively simple integrals, where Schumann derived and solved explicitly 
two infinite series utilizing Bessel functions. Schumann's sol^ion was 
for a liquid flowing through a porous solid, however, dimensional analy- 
sis may be used to make his analysis applicable to gasses [24]. Hausen's 
derivation was similar to Schumann's, and in addition Hausen also solved 
the equations for cyclic heating and cooling. Nusselt *s derivation was 
more simplified than Hausen's or Schumann's, but it is worth noting that 
Nusselt was the only one who considered the problem as part of a regenera- 
tor theory [10]. A summary of the single blow theory as it applies to 
a gas flowing through a porous solid is presented below. 

B. THEORY 

A homogeneous, porous solid has a fluid flowing through 
it. Both fluid and solid are at the same constant, uniform temperature. 
After a step change in the entering fluid temperature, the temperature of 
the solid and the fluid are to be determined as functions of time and 
position (along the flow length of the solid). 
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starting with an energy balance in the matrix: 




Assumptions : 

(1) The fluid properties are temperature independent. 

(2) The fluid flow is steady. 

(3) The porous solid is homogeneous. 

(4) The thermal conductivities of the solid and the fluid are 
infinite in the direction perpendicular to the fluid flow. 

(5) The thermal conductivity of the fluid is zero in the flow 
direction. 

Boundary conditions: 

(1) The matrix is initially at uniform temperature. 

(2) At time zero, there is a step change in the entering fluid 
tempera tu re . 

(3) The matrix boundaries are adiabatic. 
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Heat rates from the above energy diagram are: 



(1) Heat absorbed by the solid: 

(2) Heat transferred to the solid by convection: 

(3) Heat transferred from the fluid: 

(4) Heat transferred in the solid by conduction: 



h b (tp 



The net resulting energy balances for fluid and solid are: 

(1) Fluid; mCp ^^cIk + hb (t F "ts)dx - O (2-1) 

(2) Solid: PsAs Cs - koAs <2-2) 

Now introduce the following dimensionless parameters: 



(1) Time parameter 

T* iid— 

V/s 

Where: 




(2-3) 



unit conductance for convective heat transfer 



A = matrix heat transfer area (sq ft) 

VsC^ = matrix capacity (Btu/deg F) 

3 = time (hours) 

Wf = mass of fluid contained within matrix (lb m) 

• 

tn = fluid mass flow rate (Ibm/hr) 

^ = distance from matrix inlet in flow direction (ft) 

L = matrix flow length (ft) 



Introducing the fluid specific heat, and regrouping, (2-3) becomes 

HA ^ ^ KA x,_ V/r^p 






(2-4) 



The term 



V/f Cp 



1^ Cp L Ws ^6 

is very small since the thermal capacity of the matrix 



is much greater than the thermal capacity of the fluid contained within 
it. Thus, the dimensionless time parameter can be reduced 'to: 
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(2-5) 



T ~ 



hA 

V/5C5 



B 



(2) Position parameter ; 

■2 « h A X 

\fi^Cp L 

Since the Number of Transfer Units, NitU - : 



ictioi 

X- 



(2-6) 



(2-7) 



(3) Conduction parameter; 

ksAs 



Where; 

k 



s 



s 



m cpL 

matrix thermal conductivity (Btu/hr sq ft deg F/ft) 
solid matrix cross sectional area Available for thermal 



(2-8) 



conduction (sq ft) 

Introducing the above dimensionless parameters into the heat balance 
equations for the fluid and solid, and rearranging yields: 

Fluid: (2-9) 

Solid: tp-ts ■^XMU(’|^ <2-10) 

If thermal conduction in the solid is assumed zero in the direction 
of flow, the above equations simplify: 

ts-tp (2-11) 

^ (2-12) 

^ - tF-ts 

From equations (2-11) and (2-12), Schumann, utilizing Bessel func- 
tions, developed his temperature -time-distance relationships: 



^-r A 



n »1 

ts-t; . 

h-0 



JC'irzr 



(2-13) 



(2-14) 



Plots of these functions versus Z appear in Figure 1. Using these 
curves to determine Ntu from experimental data would be a tedious task 
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involving infinite series computations and curve comparisons. However, 
in his maximum slope technique, Locke [15] differentiated the theoretical 
curves, plotting slope versus "^/Ntu for various values of Ntu. Locke's 



expression for slope is: 











The downstream fluid temperature is taken at x = L so that 21 be- 
comes Ntu and t^ = Locke also evaluated the maximum slope of the 

temperature curve for each Ntu, and plotted maximum slope versus Ntu, a 
single curve. Figures 2 and 3 show Locke's cuirves. 

Both Schumann and Locke assumed thermal conductivity of the solid, 
in the flow direction, to be zero. The effect of longitudinal thermal 
conductivity is small in the high Reynolds number range, but it is con- 
siderable and should be taken into account in the low Reynolds range. 

To account for longitudinal thermal conduction, equations (2-9) and 
(2-10) must replace equations (2-11) and (2-12). Thus, the maximum 
slope will depend upon both Ntu and X . Howard [9], utilizing a 
digital computer and finite difference technique, tabulated and plotted 
Ntu versus maximum slope for various values of X . Howard's curves 
appear in Figures 4 through 6. By means of a properly designed experi- 
mental apparatus which satisfies Howard's idealizations and boundary 
conditions, the maximum slope of the heating or cooling curve and the 
conduction parameter can be determined. Ntu can then be found using 
Howard's data shown in Figure 5. 

In the vicinity of Ntu = 2, it can be seen from Figure 6 that there 
can be considerable error in Ntu for a small error in maximum slope. 
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In addition, it is more difficult to accurately determine maximum slope 
in the high Reynolds number region, as the maximum slope occurs earlier 
in time as Ntu decreases, and the determination of the maximum slope 
becomes difficult, having large uncertainties due to the inherent time 
response of the instrumentation* 

Another theory for determining the heat transfer characteristics of 
a porous solid, based upon a sinusoidally varying fluid temperature at 
the matrix inlet, was developed by Bell and Katz [2]* Starting as be- 
fore with an energy balance in the matrix: 




Assumptions: 

(1) All cross sections of the solid normal to the flow direction 
are uniform* 

(2) The fluid velocity is uniform throughout the matrix* 

(3) The fluid properties are temperature independent* 

(4) Thermal conductivity of the solid and fluid is infinite in the 
direction normal to fluid flow* 
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(5) Thermal conductivity of the fluid and solid is zero in the 



direction parallel to the fluid flow. 

( 6 ) The matrix boundaries are adiabatic. 

(7) The convective heat transfer coefficient remains constant. 
Heat rates from the above heat balance are: 

(1) Heat absorbed by the solid: 

(2) Heat transferred to the solid by convection: 

(3) Heat absorbed by the fluid; ^ 

(4) Heat transferred in the solid by conduction: o 

The net resulting energy balances for fluid and solid are: 






(1) Fluid: rnCp (2-16) 

(2) solid: PsAsCs -hb 



Inlet and initial conditions are: 

(1) Inlet: "tf (Oj&) ' £ ^ -f A *5'^ (2-18) 

(2) Initial: £5 (o) ^ t^c (2-19) 

(3) Initial: tfr(XjO)^ to (2-20) 

where ^ is the amplitude of the sinusoidally varying temperature. 

Now introduce the following dimensionless parameters: 

(1) Time ; ^ (2-21) 

(2) Distance: ^ (2-22) 

(3) Frequency: oO - (2-23) 

(4) Flow: (2-24) 

WsCs COo 

(5) Heat transfer: UA (2-25) 

Where: 



U = flow velocity 
0 = time 

L = matrix flow length 
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(ft/sec) 

(sec) 

(ft) 



X ■ distance from matrix inlet in flow direction (ft) 

angular frequency of temperature oscillation (rad/sec) 

Introducing the above dimensionless parameters into equations (2-16) through 
(2-20) and rearranging yields: 

(1) 1^ ^ Nitu(tLF-t5) (2-26) 

(^F-U ) = O (2-27) 

(3) tr(o,^) ^ tyy, ^ (2-28) 

(4) tsCo) ^ As„ (2-29) 

(5) tf (^,0) - to if) (2-30) 

This system has the periodic solution: 

izf Ae sih (co^ -jiSf) (2-31) 



t 

Where: 






5 ( f ; e M£tA Si'n (^4o"t cos (40“^ (2-32) 









/ -t- (^r1 Mta)’- 
M 

I 4 



o> 



A simple quantity to determine experimentally is R^, the ratio 
of the amplitude of outlet and inlet gas temperatures. Straightforward 
calculation shows: 



Ra- 



Ae 



- 6 



Ni. U, 

Miu)^ 



(2-33) 



A plot of Ntu versus for various values of M is shown in Figure 7. 
The same experimental apparatus can be used for this cyclic operation 
as is used for the maximum slope method. The only addition needed is 
a method to sinusoidally vary the temperature of the air entering the 
matrix. 

Thus, heat transfer data can be determined for a wide range of 



21 



Reynolds numbers utilizing two testing methods on a single experimental 



apparatus; the single-blow technique producing heat transfer data for 
low Reynolds numbers, and the cyclic technique producing data for high 
Reynolds numbers. 
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3. EXPERIMENTAL TECHNIQUES 



A. Single Blow Technique 

For Reynold's numbers up to approximately 1,000, the results 
of Howard's work [9] were used in determining the heat transfer data, 

Ntu and j, appearing in this report. In order to utilize Howard's 
theoretical data, the experimental apparatus had to conform to the ideal! 
zatlons stated in reference [9]. 

Idealizations: 

(1) The fluid flowing through the matrix remains steady and 
uniform in velocity and temperature at any cross section. The matrix 
thermal conductivity is infinite in the direction normal to fluid flow, 
and finite in the direction parallel to fluid flow. Thus the problem is 
one-dimensional in space. 

(2) The thermal capacity of the matrix is large compared with 
the fluid contained within it. This restricts the fluid to a gas, and 
means that there will be no time dependent terms in the equations for 
the fluid. 

(3) The thermal properties of the fluid and matrix are cons- 
tant and uniform. 

(4) The convective heat transfer coefficient is some suitable 
average and remains constant. 

(5) At time zero, the change in fluid temperature is a step 
change, with the matrix and its entrained fluid initially at uniform 
and equal temperatures. 

The test apparatus utilized for this report satisfies the above 
idealizations as described below: 
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(1) A square cross section contracting cone and screen type 
flow straighteners provided the necessary uniform velocity profile. In 
addition, the matrix test section was located as far upstream as possible 
to take advantage of this uniform velocity profile. The flow rate 
measuring apparatus was located downstream. The heater configuration, 

Figure 10, insured a uniform temperature profile due to even distribu- 
tion of heat over the flow cross section. The thermal conductivity in 
the direction parallel to fluid flow is taken care of by the use of the 
conduction parameter \ . 

(2) The use of a gas as the fluid flowing satisfies idealiza- 
tion (2). 

(3,4) A small temperature variation will insure uniform solid 
and fluid physical properties. At 80®F, a temperature variation of + 

10®F will result in an air viscosity variation of approximately + 1.5%. 

The changes in specific heats and thermal conductivities are negligible 
for a Hh 10®F change in temperature. 

(5) Various methods have been devised to produce a step change 
in matrix inlet temperature. The electric heater arrangement used in 
this report appears to produce the closest approximation to the desired 
step change with a minimum effect on flow properties [24]. 

A schematic diagram of the experimental apparatus is shown in 
Figure 8. Photographs of various components are shown in Figures 9 
through 14. The system is comprised of an entrance section with a contract- 
ing cone and a screen flow straightener, a heat section consisting of 28 
resistance type wire heaters made of nichrome wire in grid fashion, a 
matrix test section, an orifice metering section, a flow control valve 
and prime mover, and pressure and temperature measuring and recording 
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systems* Each of these components are further described in Appendix 
A. The experimental procedure for the single blow technique is as 
follows: 

Air is drawn through the system by means of a 30 horsepower Spencer 
Turboblower* The heater section is designed to heat the air to approxi- 
mately twenty degrees above ambient temperature* The mass flow rate is 
measured by means of an ASME orifice meter with removable orifice plates* 
Flow friction data is obtained by measuring the isothermal pressure 
drop across the matrix and the mass flow rate through the orifice meter* 
Heat transfer data is obtained as follows: the heaters are energi- 

zed to heat the air to approximately twenty degrees above ambient tempera- 
ture* When the matrix has been heated to a uniform temperature, the 
heaters are quickly de-energized and the time -temperature history of the 
air leaving the matrix is recorded on a Minneapolis Honeywell **Brown** 
strip chart recorder. 

The following data is recorded for each test run: 



(1) Atmospheric pressure Patm 

(2) Static pressure upstream of matrix 

(3) Pressure drop across the matrix A P 

m 

(4) Static pressure upstream of orifice P 

o 

(5) Pressure drop across the orifice AP 

o 

(6) Temperature upstream of orifice t^ 

(7) Orifice diameter d 

o 

(8) Time- temperature cooling curve 

(9) Recorder chart speed 
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B. Cyclic Technique: 

To obtain heat transfer data by means of the cyclic tempera- 
ture variation method, the same testing apparatus can be used, the only 
modification being in the power input to the heaters. The idealizations 
set forth by reference [2] are the same as those imposed by reference [9] 
and shown on the previous pages with two exceptions: 

(a) The thermal conductivity of the solid is zero in the 
direction parallel to the fluid flow. 

(b) The temperature of the matrix inlet varies sinusoidal- 
ly. To satisfy idealization (a), this testing method is restricted to 

the high Reynolds number ranges ( > 500) , where Ntu is small (Ntu ^ 5) , 

and the effect of longitudinal thermal conduction is small and can be 
neglected. (See Figure 15 and Chapter 5). 

To satisfy idealization (b), sinusoidally varying power is ap- 
plied to the heaters by means of a 1.5 KW Amplidyne receiving field excita- 
tion from a variable frequency function generator. This system is des- 
cribed further in Appendix B. 

To obtain heat transfer data, the temperatures of the matrix inlet 
and exit are recorded on the strip chart recorder. The amplitudes of 
the cyclic variation of these two temperature curves are then measured 
and Ntu is determined from the ratio of the amplitudes and the dimension- 
less parameter, M. Figure 15 shows a sample chart trace for a cyclic 
test. 

The same data are recorded as for the maximum slope technique. 

Data reduction relationships for both of these transient testing tech- 
niques are described in Appendix B. 
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4. 



DESCRIPTION OF MATRICES 



The surfaces investigated in this report are of two major varieties 
modified rectangular passage and triangular fin. 

The modified rectangular passage surfaces, previously examined by 
Bannon [1], and Piersall [22], are constructed of formed sheets separat- 
ed by flat splitter plates as shown below. The matrix approximates one 
made up of rectangular passages. The rectangular passages have an aspect 
ratio, oc , of approximately 7: 



One parallel plate matrix is made from .0022 inch thick perforated 
nickel and the other from .002 inch thick solid nickel. The perforations 
are of type 160/40 TV and are described further in Appendix C. 

The five triangular fin matrices, provided by Solar, A Division of 
International Harvester, San Diegp, were of varied material and method 
of construction. Geometric and physical data for these matrices appear 
in Figures 16 through 22. Two were made from nickel; one had fins made 
from .005 inch thick solid material, the other had fins made from type 
80/20 T perforated nickel material .0047 inches thick. The other three 
triangular fin matrices were made from .005 inch thick type 430 stain- 
less steel, and had no perforations. Of these stainless steel matrices, 
one had fins with very sharp bends, (all others had approximately sinusoi 
dal corrugations), one had maximum fillet braze joints where fins and 
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splitters touched, (all others were not bonded by any means) , and the 
last stainless steel surface was similar in shape and construction to 
its solid nickel finned counterpart* All of the triangular fin matrices 
had a fin height of .1 inches and 10 fins per inch. All of the matrices 
measured 3-1/16 inch high X 3-1/16 inch wide. The modified rectangular 
passage matrices were 2 inches in flow length, while the triangular 
finned ones were 3 inches long. Geometric and physical data for all 
matrices appear in Figures 16 through 22. 
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5. 



PRESENTATION OF RESULTS 



Heat transfer and flow friction data for all surfaces tested are 
shown in TABLES I through VII and Figures 23 through 35. In Figures 
23 through 31 the Colburn heat transfer modulus, j, and Fanning fric- 
tion factor, f, are plotted versus Reynolds number. The Reynolds number 
was based on hydraulic diameter. In Figure 32, all surfaces are compared 
by a plot of j/f versus Reynolds number, and Figure 33 is a plot, for 
all surfaces, of heat transfer power versus flow friction power on a 
unit area basis. 

Figure 34 is a plot of Colburn j versus Reynolds number for the brazed 
stainless steel triangular fin surface. The J values for this plot were 
obtained by means of both the single-blow technique, and the cyclic tempera- 
ture variation technique described in Chapters 2 and 3. For N < 1000 
the single-blow technique was employed, and for N > 500 the cyclic tech- 
nique was used. This was the only surface tested by means of the cycli- 
cal temperature variation technique. 

Figure 15 shows a plot of Colburn j versus Reynolds number for the 
perforated nickel rectangular passage matrix. This plot shows the ef- 
fects of varying the longitudinal conduction parameter, X . One curve 
results from assuming X = 0, the second curve results from the 
calculated in reference [22] and described in Appendix C, and the third 
curve results from an electrical analogy method calculation of X which 
is also described in Appendix C. 

Figure 35 is a plot of Colburn j versus Reynolds number for the 
perforated nickel triangular fin matrlce. the data for this curve was 
obtained by means of both the single blow technique and a steady state 
steam to air technique described in reference [29]. 
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In Figure 36, experimental and theoretical data are compared for 



the brazed triangular fin matrix. 



i 
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6 . 



EXPERIMENTAL UNCERTAINTIES 



The experimental uncertainties for the heat transfer and friction 
data appearing for this report were determined by the method described in 
reference [13] and shown in Appendix D. These uncertainties arise from 
possible errors in: 

(1) Physical constants - These physical constants were 
obtained from references [6], [7], and [8], and their uncertainties are 
estimated to be as follows: 



ks ± -5^ 

Cs + . 5% 

Cp + 

I^Pr ±2-07. 

yM + 1.0% 

(2) Geometric measurements - These errors are due to con- 
struction inaccuracies and errors in linear measurement and are esti- 
mated to be as follows: 

Ac , A , As ± 1.0% 

u + .5% 

Ws negligible 

(3) Instrumentation - Assuming adequate manufacturer's cali- 
bration a mancxneters and thermocouple wire, instrumentation errors will ^ 
be assumed to be due to instrument fluctuations (manometers) , and size 
of scale divisions (manometers and potentiometer) . 

Pressures : 

Ps + 1.0% 

APm ± 1*07. 
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Po ± 1-07- 

APo ± 1-07. 

Patm negligible 

The error for the manometers depends upon the range of the pressure 
being read. In the range 1 to 3 inches of water, the error is .17% to 
.5% in the range 3 to 5 inches of water, the error is 1.0% to 1.7%. 

Temperatures : ^ 

+ 1,5 c/«g 

Maximum Slope : 

The measurement of maximum slope involves the drawing of 
a tangent and the linear measurement of three quantities X , y , and 
in , as described in Appendix A. The percent error in maximum slope 
is estimated as + 2.0%. Entering Figure (7) with this value at = 3.0 
and ^ = 0, obtain Ntu error as + 7.2%. At Ntu = 25 and X = .05, Ntu 
error is + 10.0% for a 2.0% error in maximum slope. 

A summary of extreme values of uncertainty intervals calculated by 
the method shown in Appendix D is given below: 

± 1 . 0 % 

Nr + 2-37. 

j ± W’l 

f ± 
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7. 



DISCUSSION OF RESULTS 



Heat transfer and flow friction results appear in Figures 23 through 
36, with comparisons of all surfaces in Figures 31 through 33. 

In general, the triangular fin surfaces all had similar j and f 
curves; the nickel matrices having higher j and j/f values than the stain- 
less steel matrices. The high j and j/f for the nickel surfaces are due 
chiefly to the high thermal conductivity of nickel - 36.0 compared to 
12.8 for 430 stainless steel (at 80®F). The j,f, and j/f characteristics 
for the triangular fin surfaces are in agreement with the theory presented 
in reference [16] for laminar flow through idealized cross sections (see 
Figure 36). 

The three stainless steel cores were quite similar in cross section, 
however, one might expect to find small differences in their heat trans- 
fer and flow friction characteristics due to slight differences in geo- 
metry. The rounded comers produced by the sine wave corrugations in 
two of the surfaces should reduce the effective surface area for heat 
transfer and wall friction as mentioned in reference [16]. The experi- 
mental results seem to agree in part with the above - the sharp cornered 
triangular passage geometry yielded higher friction factors than the 
sinusoidal fin geometry, with the brazed sinusoidal fin core producing 
intermediate friction factors. The heat transfer characfieristics, how- 
ever, showed the j curve for the brazed core slightly higher than the j 
curves for the unbrazed ones. There are two reasons for this: First of 

all, the heat transfer area calculations for the brazed core took into 
account the loss of metal surface area due to brazing. In the unbrazed 
sinusoidal fin core, a portion of surface area similar to the **lost area” 
was relatively ineffective for heat transfer, but not accounted for in 
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calculating j* Secondly, the fin shape and matrix construction for the 
brazed core actually produced a geometry that is a cross between a tri* 
angular passage and a rectangular passage, with a rectangular passage 
having higher values of J and j/f than a triangular passage as shown in 
reference [16]. 

It should be noted that these three stainless steel matrices yielded 
heat transfer and flow friction results that differed only slightly, in- 
dicating that small deviations in manufacturing procedure may have little 
effect upon heat transfer and flow friction performance. 

The rectangular passage cores had heat transfer and friction charact- 
eristics similar to the triangular fin cores; however, their heat trans- 
fer curves (j versus N ) were less steep, the perforated material pro- 
ducing a markedly higher curve than the solid material. The j/f flow 

area "goodness” factors for these surfaces increased as N increased, 

R 

while j/f for the triangular fin surfaces decreased with increasing 
Reynolds number (See Figure 32). On the h versus E plot. Figure 
33, rectangular passage surfaces showed definite superiority over the 
triangular fin surfaces, with the perforated material producing the best 
results. 

Figure 15 demonstrates the influence of conduction parameter upon 
the heat transfer modulus, j. Curves such as those shown dotted can 
result from an improper assumption ( X = 0 in low Reynolds range) , ap- 
proximations for X (Appendix C) , or an inaccurate value for k . 

s 

Referring to Figure 7, it can be seen that there are areas where a 
small error in Amplitude Ratio, R^, would result in a large correspond- 
ing error in Ntu. However, the frequency of temperature oscillation, 

> can be chosen so as to produce a flow parameter, M, (M = me 
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<■ • that will result in Ntu being determined with minimum error for 

a given error in R^. For instance, at Ntu = 1.5, should be chosen so 
that M is less than .35, since values of M near .7 in this case result 
in a high uncertainty of Ntu. This flexibility is not inherent in the 
single-blow technique (where Ntu « 2). 

Figure 35 shows heat transfer data for the brazed triangular fin 
matrix. Both the single-blow and cyclic techniques were employed, thus 
producing data up to a Reynolds number of 1850 with an area of overlap 
between Reynoldd numbers 380 and 560. There seems to be fairly good 
correlation between the two testing methods with the exception of the 
data points at Reynolds numbers 450 and 560. These points were the re- 
sult of employing the single-blow technique in the area of high N and 
low Ntu (1«65 and 1.44), where Ntu is subject to error as explained in 
Chapter 2. The data from the cyclic technique was also subject to some 
error due chiefly to the cyclic power input to the heaters. As can be 
seen on the sample chart trace in Figure 37, there was some distortion 
in the temperature waveforms making them slightly different than true 
sinusoids. This distortion is due to inaccurate signal reproduction in 
the amplidyne. There was also a problem in instnomentation - the tempera- 
ture waveforms were recorded separately by means of the same instrument 
instead of simultaneously. Bell and Katz [2] utilized a dual channel 
instrument recording both temperatures simultaneously, reversing the in- 
put connections midway through each run so as to record each temperature 
on both channels, thus eliminating the effects of unequal channel sens!** 
tivities. The thermocouple configuration described in Appendix B does 
not allow for simultaneous temperature measurement and recording. 
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In Figure 35 data from the single-blow technique is plotted with 
data from a steady state stream to air technique. In both tests the 
identical type 80/20T perforated nickel triangular fin surface was in- 
vestigated. The steam to air test was conducted on a cross flow heat 
exchanger with the surface to be tested on the air side. The core con- 
figuration for the single-blow test is described in Chapter 4 and Figure 
17. Theoclitus in reference [28] utilized a dual testing technique similar 
to this in determining fin efficiency, , and total surface temperature 
effectiveness 7^^ , for various surfaces. Theoclitus found for simple 
finned surfaces, that the theories for fin efficiency and total surface 
temperature effectiveness are reliable since he obtained a close corre- 
lation of data from the two methods of testing. For the steady state 
technique, theoretical values of ^^and are used in determining the 
heat transfer characteristics, Ntu and j, while and ^ have values 
of unity when the transient technique is employed, due to lack of a 
temperature gradient in the fin in the direction perpendicular to fluid 
flow (See Chapter 3) . 

Figure 35 indicates a fairly good correlation of data for the two 
testing methods. 

Conclusions : 

(1) For matrices of the same geometry and material, small devia- 
tions in manufacturing procedure have little effect upon overall heat 
transfer and friction performance. 

(2) The single-blow and cyclic transient testing techniques 
utilized for this report are reliable in that the experimental results 
obtained agree with results obtained by means of steady state testing 
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and with theory. 

(3) The 
in producing heat 



single-blow and cyclic techniques supplement each other 
transfer data over a wide range of Reynolds numbers. 



8, RECOMMENDATIONS FOR FURTHER WORK 



Refinement of the cyclic method of testing is needed both in theory 
and in testing apparatus. The theory could be extended to include the 
effect of longitudinal thermal conduction, this involving the use of 
another non-dimensional parameter in place of M, and a theoretical solu- 
tion using methods similar to Howard's [9] or Moreland's [20]. 

There are a number of areas for improvement in the test apparatus: 

(1) A modification is needed if higher flow rates are to be 
employed. Either a larger pipe diameter at the orifice 
meter, or another means of measuring flow rate is needed 
to increase the flow rates employed. 

(2) A more reliable power source is needed for the heaters so 
that a truer sinusoidal temperature variation results. 

(3) The temperature measuring system could be modified so 

that the matrix inlet and exit temperatures can be recorded 
simultaneously. 

It is also recommended that more surfaces be made available for 
testing by both transient and steady state methods so that data is avail- 
able for a wider range of Reynolds numbers, and the reliability of the 
two methods may be compared. 



38 



9. 



BIBLIOGRAPHY 



(1) Bannon, John M. "An Experimental Determination of Heat Transfer and 

Flow Friction Characteristics of Perforated Material for Com- 
pact Heat Exchanger Surfaces." Master's Degree Thesis U. S. 

Naval Postgraduate School, Monterey, California. 1964. 

(2) Bell, J. C. and Katz, E. F. "A Method for Measuring Surface Heat 

Transfer Using Cyclic Temperature Variations", Heat Transfer 
and Fluid Mechanics Institute . Berkeley, California, 1949. 

(3) Carslaw, H. S. and Jaeger, J. C. Conduction of Heat in Solids , Second 

Edition, Oxford; Clarendon Press, 1959. 

(4) Coppage, J. E. "Heat Transfer and Flow Friction Characteristics of 

Porous Media", T.R. No. 16, Department of Mechanical Engineering , 
Stanford University . Stanford, California, December 1, 1952. 

(5) Creswick, F. A. "A Digital Computer Solution of the Equations for 

the Transient Heating of a Porous Solid Including the Effects of 
Longitudinal Conduction", Industrial Mathematics , Volume 8, 

1957. 

(6) Elridge, E. A. and Deem, H. W. Report on Physical Properties of 

Metals and Alloys from Cryogenic to Elevated Temperatures . 
Philadelphia: American Society for Testing Materials, 1961. 

(7) Goldsmith, A., et. al . Handbook of Thermophysical Properties of 

Solid Materials . Revised Edition, Volume I, Elements, New 
York: The MacMillan Company, 1961. 

(8) Hllsenrath, J. , et al . Tables of Thermodynamic and Transport Properties 

of Air, Argon. Carbon Monoxide, Hydrogen, Nitrogen, Oxygen, and 
Steam . New York; Permagon Press, 1960. 

(9) Howard, C. P. "The Single Blow Problem Including the Effects of 

Longitudinal Conduction", ASME paper . Number 64-GTP-ll 

(10) Jak^^' M. Heat Transfer , Volume II, New York: John Wiley and SonS, 

Inc., 1957. 

(11) Kays, W. M. "Loss Coefficients for Abrupt Changes in Flow Cross 

Section with Low Reynolds Number Flow in Single and Multiple 
Tube Systems", TR Number 9, Department of Mechanical Engineering , 
Stanford University . Stanford, California, 1950. 

(12) Kays, W. M. and London, A. L. Compact Heat Exchangers . Second Edition. 

New York; McGraw-Hill Book Company, Inc., 1964. 

(13) Kline, S. J. and McClintock, F. A. "Describing Uncertainties in 

Single-Sample Experiments," Mechanical Engineering , January 1953, 
pp ; 3-8. 



39 



(14) Lambertson, T. J. "Performance Factors of a Periodic Flow Heat 

EJtChanger", Trans. ASMS , Volume 80, 1958. 

(15) Locke, G. L. "Heat Transfer and Flow Friction Characteristics of 

Porous Solids", TR. No. 10, Department of Mechanical Engineer - 
ing, Stanford University , Stanford, California, June 1, 1950. 

(16) London, A. L. "New Developments in Compact Heat Exchangers -Design 

Theory, Surfaces, and Applications", lecture presented at the 
Sixth National ASME-AIC E Heat Transfer Conference , Boston, 
Massachusetts, August, 1963. 

(17) London, A. L. , Biancardi, F. R. , and Mitchell, "The Transient Response 

of Gas Turbine Plant Heat Exchangers, Regenerators, Intercoders, 
Precooters and Ducting, Trans . ASME , Volume 81, pp. 43, 1959. 

(18) McAdams ) Heat Transmission , Third Edition, New York: McGraw-Hill 

Book Company, Inc., 1954. 

(19) Mondt, J. R. "Effects of Longitudinal Thermal Conduction in the 

Solid on Apparent Convection Behavior, with Data for Plate Fin 
Surfaces", 1961 International Heat Transfer Conference , 

Boulder, Colorado, Paper No. 73, Proceedings ASME. 

(20) Moreland, F. E. "Solution of the Single Blow Problem, with Longi- 

tudinal Conductivity by Numerical Inversion of LaPlace Trans- 
forms", Master's Degree Thesis, U. S. Naval Postgraduate School , 
Monterey, California, 1964. 

(21) Murdock, J. W. "Tables for the Interpolation and Extrapolation of 

ASME Coefficients for Square-Edged Concentric Orifices", ASME 
paper. Number 64-WA/FM-6. 

(22) Piersall, C. H. "Experimental Evaluation of Several High Performance 

Surfaces for Compact Heat Exchangers". Master's Degree Thesis, 

U. S. Naval Postgraduate School, Monterey, California, 1965. 

(23) Power Test Code Supplements, Instruments and Apparutus , PTC-19.5.4 

Flow Measurements, Chapter 4, New York: American Society of 

Mechanical Engineers, 1959. 

(24) Pucci, ?. F., Howard, C. P. , and Piersall, C. H. , "The Single-Blow 

Transient Testing Technique for Compact Heat Exchanger Sur- 
faces," ASME Paper , Number 66-GT-93. 

(25) Raven, F. H. "Automatic Control Engineering^^ New York: McGraw-Hill 

Book Company, Inc., 1961. 

(26) Schumann, T. E. W. "Heat Transfer: A Liquid Flowing through a Porous 

Prism", Journal of the Franklin Institute, Volume 208, 1929, 
pp 405-416. 



40 



(27) Solar Aircraft Company, Engineering Report ER 1221, Regnerator core 

Module Tests-T-600 Engine, 18 April, 1962. Rev. A; 25 October, 
1962, AD 290234. 

(28) Theoclltus, G. and Ecbrlch, T. L. "An Experimental Technique for 

Determining the Effectiveness of Extended Surfaces", presented 
at the Seventh National Heat Transfer Conference . Cleveland, 
Ohio, August, 1964. 

(29) Ward, J. P. "Steady State Steam to Air Testing Facility for Compact 

Heat Exchangers". Master's Degree Thesis, U. S. Naval Postgrad - 
ate School, Monterey, California, 1965. 



41 



bat* 




TBMPERATOHE IS A PWCTION OP T KSD i OP A PLTJID 
PLOimO THROUGH AK INITIALLY UNIFORM TEMPERATURE 
POROUS SOLID 



»» |4* 



X.O 



• 9 

• 6 
.7 

•4 
• 3 
.Z 

a 

0 




SOLID TEMPERATURE AS A FUNCTION OP T AND i 
OF A P0R003 SOLID IKITIALLX AX A TOIPDIUI 
TKMPERATUiE AMD THROOOH WHICH A FLUID FLOWS 



Figure 1. Schumann's Curves 
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Figure 3. Locke's Curves 
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Ntu versus Maximum Slope 
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Figure 5. Ntu versus Maximum Slope 
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Figure 6. Ntu error versus Ntu 
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Figure 7. Ntu versus Amplitude Ratio 
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Figure 15. Effects of Longitudinal Conduction Parameter 
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Matrix Material 
Specific Heat, c ^ 

Thermal Conductivity, k ' 
Material Thickness 
Total Heat Transfer Area 
Frontal Area, 

Conductive Area, A 

s 

Free Flow Area, A 
Matrix Flow Length, L 
Matrix Volume, V ™ 
Matrix Density, 
Compactness, 

Porosity, p 
Hydraulic Diameter 



SOLAR #1 



(Btu/lbm deg F) 
(Btu/hr ft deg F) 
(in) 

,A (sq ft) 

(sq ft). 

(sq ft) 

(sq 'ft) 

' . (ft) 

(cia ft) 

(Ibm/cu ft) 

(sq ,ft/cu ft) 

(fb) 



Nickel 

.106 

36.0 

.005 

12.768 

.07022 

.01078 

.05944 

.'24917 

.0175 

86.1378 

729.7026 

.8464 

.0046395 




» 

Figure 16 Geometric and Physical Propertied 
for Solar #1 
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Matrix Material 


SOLAR #2 


Nickel 


Specific Heat, c 


(Btu/lbm deg F) 


.106 


Thermal Conductivity, k 


(Btu/hr ft deg F) 


36.0 


Material Thickness 


(in) 


.005 


Total Heat Transfer Area 


,A (sq ft) 


12.093 


Frontal Area, A^^ 
Conductive Area, A 


(sq ft) 


.07022 


(sq ft) 


.01034 


Free Flow Area, A * 


(sq ft) 


.05987 


Matrix Flow Length, L 


(ft) 


.25583 


Matrix Volume, V “ 


(cu ft) 


.01797 


Matrix Density, 


(Ibm/cu ft) 


70.3099 


Compactness, 


(sq ft/cu ft) 


673.143 


Porosity, p 


* 


.85281 


Hydraulic Diameter 


(ft) 


.0046272 


Perforated Fin Material 

■ / ' 


.0047" thick 


80/20T 



Figure 17 Geometric and Physical Properties 
for Solar #2 
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SOLAR #3 


• 


Matrix Material 




430 Stainless Steel 


Specific Heat, c 


(Btu/lbm deg F) 


.11 , 


Thermal Conductivity, k 


(Btu/hr ft deg F) 


12.8 


Material Thickness ® 


(in) 


.005 


Total Heat Transfer Area 


•,A (sq ft) 


14.3588 


Frontal- Area, A^^ 1 

Conductive Area,, A • 


(sq ft) 


,.06803 


(sq ft) 


.01325 


Free Flow Area, A 


(sq ft) 


i 05478 


Matrix Flow Length, L 
Matrix Volume, V ™ 


(ft) 


>24833 ' 


(cu ft) 


.01689 


Matrix Density, 


(Ibm/cu ft) 


94.2656 


Compactness^ 


(sq ft/cu ft) 


849.8962 


Porosity, p 
Hydraulic Diameter 


(ft) 


.80517 ' , '■ 

.0037895 



Figure 18. Geometric and Physical Properties 

For Solar #3 




. SOLAR #4 



Matrix Material 
Specific Heats c 
» Thermal Cdiaduetivitys k 
Material Thickness ® 
Total Heat Transfer Area 
Frontal Area, A. 
ConductiTe Area, A 
Free Flow Area’, A ® 
Matrix Flow LengtH, L 
Matrix Volume, V 
Matrix Density, ™ 
Compactness, ™ 
Hydraulic Diameter 



430 stainless steel 
(Btu/lbm degF) .11 

(itu/hr ft deg F) 12.8 



( in) 

,A (sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(Ibm/cu- ft) 

(sq ft/cu ft) •• 
■ (ft) 



.005 

.14.27694 

.06713 

.01229 

.05484 

.24417 

.01638 

88.5918 

.81690 

.00375168 



V Figure 19. 



Geometric and Physical Properties 
for Solar #4 



■i 



i 



60 








SOLAR #5 



Matrix Material 
Specific- Heat, c 
Thermal Conductivity, k 
Material Thickness 
Total Heat Transfer Area ,A 
'Frontal Area, A. 

If 

Conductive Area, A 
Free Flow Area, A ® ' 

Matrix Flow LengtH, L 
Matrix Volume, V “ 
Matrix Density, ™ 
Compactness, 

Porosity, p 
Hydraulic Diameter 



(Btu/lbm deg F) 
(Btu/hr ft deg F) • 
(in) 

(sq ft) 

(sq ft) 

(sq ft) 

(sq ft) 

(ft) 

(cu ft) 

(Ibfh/cu ft) 

(sq ft/cu ft) 

(ft) 



430 stainless steel 
. .11 
12.8 
.005 
14.4234 
.06781* 

■ .01148 
.05633 
.26417 
.0179 
81.9345 
805.13313 
' .83066 

.0041269 






Figure 20. Geometric and Physical Properties 

for Solar #5 ^ 
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MODIFIED . RECTANGULAR PASSAGE 



Matrix Material 

Specific Heats e , 

, Thermal Conductivity, k 

Material Thickness 

Total Heat Transfer Area ,A 

Frontal Area, A^^ 

Conductive Area, A^ 

Free Flew Area, A 

Matrix Flew Lengtn, L 

Matrix Volume , V . ™ 

# in 

Matrix Density 

Compactness, 

Porosity, p 

Hydraulic Diameter 



(Btu/lbm deg F) 
(Btu/hr ft deg F) 
(in) 

(sq ft) 

(sq' ft) 

(s'q ft) 

(sq ft) 

(ft) 

(cu ft) 

(Ibm/cu ft) 

(sq ft/cu ft) 

V 

(ft) 



160/40 Perforated nickel 
.106 
36.0 
.0022 
18.593 
.06953 
.00424 
.0584 
.16098 
.01159 
63.6 
1604.5 
.840 
.001929 



•f • -V. • ■ c. 

Figure 21. Geometric and Physical Properties 

for 160/40 Tf 

Modified Rectangular 
Passage 
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8* THEL.S,STARRETTCO. 



MODIFIED 


RECTANGULAR PASSAGE 




Matrix Material 
Specific Heat, c 
Thermal Conductivity, k 
Material Thickness 


(Btu/lbm deg F) 


Solid Nickel 
.106 


(Btu/hr ft deg F) 


36.0 


(in) 


.002 


Total Heat Transfer Area, A 


(sq ft) 


20.1875 


Frontal' Area, A^^ 
Conductive Area, A 


(sq ft) 


.06953 


(sq ft) 


• ' .01014 


Free Flow Area, A 


(sq ft) 


.05939 


Matrix Flow Length, L 


(ft) 


.16098 


Matrix Volume, V 


(cu ft) 


.011588 


Matrix Density *1? 


(Ibm/cu ft) 


77.7 


Compactness, 


(sq ft/cu ft) 


1742.1 . 


Porosity, p 
Hydraulic . Daimeter 


(ft) 


;854' 

.001961 



Figure 22. Geometric and Physical Properties 
for Solid Nicktj Modified 
Rectangular Passage 
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Figure 23. Heat Transfer and Flow Friction 
Characteristics Solar #1 
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Figure 24, Heat Transfer and Flow Friction 
Characteristics Solar #2 
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Figure 25. Heat Transfer and Flow Friction 
Characteristics Solar #3 
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Figure 26. Heat Transfer and Flow ‘Friction 
Characteristics Solar #4 
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Figure 27. Heat Transfer and Flow Friction 
.. Characteristics Solar #5 
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Figure 28. Heat Transfer and Flow Friction 
Characteristics 160/40 TV 
Modified Rectangular 
Passage 
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Figure 29. Heat Transfer and Flow Friction 

Characteristics Solid Nickel 
Modified Rectangular Passage 
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Figure 30. Comparison of Heat Transfer 

Characteristics for all Surfaces 
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Figure 31. Comparison of Flow Friction 

Characteristics for All Surfaces 



72 




100 



1000 



2000 



N. 



Figure 32. Comparison of Flow Area Goodness 
Factors for All Surfaces 
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Figure 33. 



Heat Transfer Power versus 
Flow Friction Power per Unit Area 
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Figure 34. Heat Transfer Data by* Means of 
Cyclic Technique Solar #3 
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Figure 35. Heat Transfer Data by Me^ns of Single Blow 
and Steady State Steam to Air Test Solar t.2 



76 




j/. 



.1 



.01 



!TTOLlM!!l!ijlll!i 
' " ^pij' 

lilili; 




HEAT TRANSFER and FLQV FRICTION DATA 
SOLAR #3 BRAZED TRIANGULAR FIN 

EXPERIMENTAL DATA, SINGLE BLOW and CYCLIC TECHNIQUES 



THEORETICAL DATA, LAMINAR FLOW, TRIANGULAR PASSAGE 



,001 



5000 



N, 



Figure 36. Comparison of Experimentally Determined 
Heat Transfer and Friction Dat^ with 
Theory Solar #3 
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Figure 37. Chart Trace From Cyclic Test 
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Figure 38. Conduction Parameter Analog 



79 





TABLE I. Summary of Heat Transfer and Friction Results 



SCLAR NC. 1 SINUSOIDAL FINS SOLID NICKEL 5 MLS 



FRICTION DATA HEAT TRANSFER DATA 





CONC 


MAX 






RUN NR F 


PARAM 


SLOPE 


NTU 


NR J 


1 35.66 .58553 


.32719 


.92760 


25.57 


35.14 .09482 


2 69. SI .3013U 


. 16643 


.93616 


23.81 


68.89 .08829 


3 1C6.C0 .20595 


. 10970 


.92300 


14.74 


104.47 .05464 


4 136.80 .15680 


.08432 


.88946 


11.21 


134.83 .04155 


5 170.53 . 12574 


.06764 


.86085 


9.38 


168.08 .03475 


6 267.73 .08029 


.04312 


.72688 


5.48 


263.87 .02033 


7 394.40 .05971 


.02925 


.63549 


3.78 


388.72 .01401 


8 523.69 .04418 


.022C3 


.56660 


2.59 


516.15 .00959 


9 642.29 .03663 


.01795 


.54626 


2.09 


633.05 .00775 


10 912.84 .02736 ■ 


.01263 


.49017 


1.42 


899.71 .00525 


11 1164.10 .02248 


.00990 


.39251 


1.06 


1147.36 .00394 


12 1293.89 .02064 


.00892 


.39769 


1 .08 


1275.26 .00400 


13 1432.01 .01873 


.008C5 


.37934 


1.03 


1411.41 .00380 


14 1557.13 .01773 


.00741 


.37751 


1.02 


1534.71 .00378 


15 1678.22 .01681 


.00688 


.30071 


.00 


1654.05 .00000 


FRCNTAL AREA 


AFR 


• 


07022 




FLCW LENGTh 


L 


• 


24917 




FATRIX DENSITY 


LB/FT3 


86. 


13783 




FREE FLOW AREA 


AC 


• 


05944 


— 


SOLID FRCNTAL AREA 


AS 


• 


01078 




POROSITY 


POR ■ 


• 


84640 




HEAT TRANSFER AREA 


A 


12 . 


76800 




HYDRAULIC RADIUS 


RH 


• 


00116 


— 


AREA DENSITY (A/VM) 


( BETA) 


729. 


70260 




CCNDUCTION PARAMETER 


RATIO (RS/RK) 1. 


OCCCO 




THERMAL CONDUCTIVITY 


KS 




36.00 




MATRIX HEAT CAPACITY 


CS 




.1060 




entrance press. loss 


CCEF. KC 


• 


93767 




EXIT PRESSURE LOSS 


COEF. KE 




70431 
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TABLE II 



Summary of Heat Transfer and FtlctiiOn Results 



SOLAR NO. 2 80/20T SINUSOIDAL FINS SOLID NICKEL SPLITTERS 5 ^«ILS 





FRICTION DATA 




HEAT 


TRANSFER DATA 








CCNC 


MAX 








RUN 


NR F 


PARAM 


SLOPE 


NTU 


NR 


J 


1 


38.67 ,53389 


. 14025 


.95455 


25.60 


38.12 


.10087 


2 


75.32 .26494 


.C72C8 1 


.03187 


21.88 


74.24 


.08623 


5 


110.81 .18811 


.04926 


.96874 


13.84 


109.21 


.05455 


4 


139. OC .15304 


.03968 


.91656 


1C. 98 


136.98 


.04330 


5 


212.84 .10594 


.02560 


.77382 


6.60 


209.78 


.02602 


6 


276.05 .08185 


.01996 


.70906 


5.24 


272.04 


.02066 


7 


421.28 .05593 


.01290 


.58821 


3.09 


415.24 


.01217 


8 


497.17 .04840 


.01096 


.54974 


2.35 


490.02 


.00926 


9 


570.53 .04258 


.00963 


.53125 


1.72 


562.26 


.00679 


10 


705.01 .03521 


.00779 


.50567 


1.52 


694.80 


.00598 


1 1 


823.83 .03072 


.00659 


.46394 


1.30 


812.02 


.00513 


12 


975.36 .02720 “ 


.00562 


.39936 


1.09 


961 .25 


.00428 


13 


1096.74 .02505 


.00500 


.38556 


1.05 


io8o.ee 


.00412 


14 


1368.82 .02142 


.00399 


.00000 


.00 


1349.10 


.00000 


15 


1538.31 .02057 


.00355 


.00000 


.00 


1516.14 


.00000 


16 


2317.65 .01848 


.00226 


.00000 


.00 


2284.25 


•00000 


FRCNTAL AREA 


AFR 


• 


07C22 






FLOW 


LENGTH 


L 


• 


25583, 






PATRIX DENSITY 


LB/FT2 


70. 


30855 






FREE 


FLOk« AREA 


AC 


• 


05989 






SCLIO FRONTAL AREA 


AS 


• 


01C34 . 






PCROSITY 


POR 


• 


85281 






HEAT 


TRANSFER AREA 


A 


12. 


093C2 






HYDRAULIC RADIUS 


RH 


• 


00116 






AREA 


DENSITY (A/VM) 


( BETA ) 


673. 


13967 






CCNDUCTIGN PARAMETER 


RATIO IRS/RK 


) 


46730 




• 


THERMAL CCNCUCTIVITY 


KS 




36.00 




- 


MATRIX HEAT CAPACITY 


CS 




. 1060 






ENTRANCE PRESS. LCSS 


CCEF. KC 


• 


07222 






EXIT 


PRESSURE LCSS 


CCEF. KE 


• 


02166 
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}. ' I TABLE III. Summary of Heat Transfer and Friction Results 



SCLAR NO. 2 SINUSOIDAL FINS SOLID 430 STAINLESS FAX ERAZE 5 MILS 



FRICTION DATA HEAT TRANSFER DATA 





CONC 


MAX 








RUN NR F 


PARAM 


SLOPE 


NTU 


NR 


J 


1 60.19 .20320 


.07438 1 


.03852 


23.30 


59.32 


.07076 


2 91.18 .13761 


.049C7 


.97702 


14.24 


89.87 


.04325 


3 121.78 . 10327 


.03673 


.891 37 


9.98 


120.03 


.03031 


4 172.54 .07288 


.02593 


.77241 


6.57 


170.06 


.01995 


5 221 .36 .06225 


.02021 


.68751 


4.82 


218.18 


.01464 


6 342.40 .04132 


.013C8 


.57362 


2.81 


337.48 


.00853 


7 456.46 .03125 


.00980 


.52398 


1.65 


449.90 


.00502 


8 567.42 .02574 


.00789 


.49362 


1.44 


559.26 


.00438 


9 672.09 .02200 


.006f6 


.45011 


1 .25 


662.44 


.00379 


10 902.44 .01737 


.00496 


.00000 


.00 


889.47 


.00000 


11 1194.66 .01378 


.00375 


.00000 


.00 


1177.47 


.00000 


12 1436.51 .01193 


.00312 


.00000 


.00 


1415.83 


.00000 


13 1921.60 .00975 


.00233 


.00000 


.00 


1893.93 


•OOOCO 


FRONTAL AREA 


AFR 


• 


06803 






FLOW LENGTH 


L 


• 


24833 






MATRIX DENSITY 


LB/FT3 


94. 


26563 




- 


FREE FLOV» AREA 


AC 


• 


05478 






SOLID FRONTAL AREA 


AS 


• 


01325 






POROSITY 


PCR 


• 


8C517 






HEAT TRANSFER AREA 


A 


14. 


35884 




- 


HYDRAULIC RADIUS 


RH 


• 


0C095 






AREA DENSITY (A/VM) 


( BETA) 


849. 


89621 






CONDUCTION PARAMETER RATIO (RS/RK 


) • 1. 


OCOOO 






THERMAL CCNDUCTIVITY’ ■ 


KS 




12.80 






MATRIX HEAT CAPACITY 


cs . 




.1100 






ENTRANCE PRESS. LOSS CCEF 


. KC 


• 


97439 






EXIT PRESSURE LOSS CCEF 


UJ 

• 




65449 
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TABLE IV. Summary of Heat Transfer and Friction Results 






SOLAR NO.M SINUSOIDAL 


FINS SOLID 


430 STAINLESS 


MAIL. 5 MILS 




FRICTION DATA 




HEAT 


TRANSFER DATA 








CONC 


MAX 








RUN 


NR F 


PARAM 


SLOPE 


NTU 


NR 


J 


1 


30.98 .37693 


. 13445 


1 .16696 


2.00 


30.53 


.00611 


2 


59.13 .20381 


.C7C44 


1.07385 


27.52 


58.28 


.08412 


Zl 


88.76 . 14272 


-.04692 


.97752 


14.10 


87.49 


.04310 


4 


1 12.87 . 11402 


.03689 


.89210 


10.01 


111.25 


.03061 


5 


171.76 .07412 


.02425 


.74849 


6.04 


169.30 


.01645 


6 


215.99 .05929 


.01929 


.68279 


4.74 


212.89 


.01449 


7 


338.77 .03983 


.01220 


.5641 1 


2.64 


333.91 


.008C6 


8 


447.77 .03100 


.00929 


.50780 


1.53 


441.35 


.00467 


9 


557.45 .02546 


.00747 


.48585 


1 .40 


549.45 


.00429 


1C 


664.38 .02172 


.00627 


.41536 


1.13 


654.85 


.00346 


n 


888.19 .01741 


.00469 


.00000 


.GO 


875.45 


.00000 


12 


1162.42 .01426 ' 


.00358 


.00000 


.00 


1145.75 


.00000 


13 


1415.34 .01239 


.00296 


.00000 


.00 


1394.95 


.OOOCO 


14 


1908.14 .01161 


.00219 


.00000 


.00 


1880.66 


.00000 


FRONTAL AREA 


AFR 


• 


06713 






FLOW 


LENGTH 


L 


• 


24417 




1 


MATRIX DENSITY 


LB/FT2 


• 

CO 

00 


59182 






FREE 


FLOW AREA 


AC 


• 


05484 






SOLID FRONTAL AREA 


AS 


• 


01229 






POROSITY 


PCR 


• 


81690 






HEAT 


TRANSFER AREA 


A 


14. 


27694 






HYDRAULIC RADIUS 


RH 


• 


00C94 






AREA 


DENSITY (A/VM) 


( BETA) 


870. 


99223 






CONDUCTION PARAMETER 


RATIO (RS/RK) 1. 


OCOOO 


- 




THERMAL CONDUCTIVITY 


KS 




12.80 






MATRIX HEAT CAPACITY 


CS 




.1100 






ENTRANCE PRESS. LOSS 


CCEF. KC 


• 


96374 






EXIT 


PRESSURE LOSS 


CCEF. KE 




66901 




— 



83 



TABLE V. 



Summary of Heat Transfer and Friction Results 



SCLAR NC.5 triangular FINS SOLID 430 STAINLESS FATL. 5 MILS 



FRICTION DATA HEAT TRANSFER DATA 







CONC 


MAX 








RUN 


NR F 


PARAM 


SLOPE 


NTU 


NR 


J 


1 


63.05 .22228 


.061 18 1 


.04806 


21.15 


62.15 


.06573 


2 


96.54 . 14597 


.03995 


.95656 


12.62 


95.15 


.03924 


3 


121 .54 . 1 1431 


.03172 


.88205 


9.54 


119.79 


.02967 


4 


183.46 .08385 


.02101 


.73707 


5.81 


180.83 


.01805 


5 


231.11 .06594 


.01667 


.66779 


4.48 


227.80 


.01393 


6 


360.72 .04488 


.01068 


.56045 


2.57 


355.54 


.00800 


7 


482.55 .03439 


.00798 


.50007 


1.48 


475.63 


.00460 


8 


600.21 .02846 


’ ■ .00642 


.47897 


1.37 


591.59 


.00426 


9 


708.74 .02472 


.00544 


.45696 


1.28 


698.57 


.00396 


10 


964.34 .01907 


.004C1 


.00000 


.00 


950.46 


.00000 


1 1 


1247.65 .01548 


.00310 


.00000 


.00 


1229.68 


.00000 


12 


1548.85 .01341 


.00250 


.00000 


.00 


1526.55 


.00000 


13 


2040.53 .01207 


.00190 


.00000 


.00 


2011.13 


•ooooo 


FRONTAL AREA 


AFR 


• 


06781 






FLOW 


LENGTH 


L 


• 


26417 






MATRIX DENSITY 


LB/FT2 


81. 


93452 




- 


FREE 


FLOW AREA 


AC 


• 


05633 






SOLID FRONTAL AREA ' 


AS 


• 


01148 






PCROSITY 


PCR 


• 


82066 






HEAT 


TRANSFER AREA 


A 


14. 


42340 


■ - 





HYDRAULIC RADIUS 


RH 


• 


0C1C3 






AREA 


DENSITY (A/VM!) 


(BETA) 


805. 


13313 






CCNDUCTION PARAMETER 


PATIO (RS/RK 


) 1. 


OCCOO 







12.80 
. nco 

.95140 

.68569 



THERMAL CONDUCTIVITY KS 
MATRIX HEAT CAPACITY CS 
ENTRANCE PRESS. LOSS CCEF. KC 
EXII_PRESSURE LOSS CCEF. KE 
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TABLE VI* Summary of Heat Transfer and Friction Results 



160/40FV P6RF0RAI6D NICKEL FINS AND SPLITTERS 





_ friction 


DAf A 




HEAT TRANSFER 


DATA __ 








CCND 


MAX 




run 


NR 


F“ 


parah 


SLOPE NiU' 


~ NR J 



1 


812.99 .02545 


.00318 


■ .55590 


2.54" 


801.26 .00635 


2 


653.01 .03104 


.00395 


.59357 


3.27 


643.60 .00816 


3 


527.22 .03751 


. 00490 


.65507 


4.36 


519.61 .01091 


” ' 4 


406.11 .04445 


.00635 


“ .69800"' 


5.15 


'400 .27". 01287 


5 


310.00 .06028 


.00832 


. 77426 


6.63 


305.54 .01650 


6 


206.30 .00782 


.01253 


. 89948 


9.62 


203.40 .02405 


' - 7 


124.11 .13378 


. 02000 


1 .06826 


15.73 ' 


122.32 .03932 


8 


69.91 .23837 


. 03694 


1.17903 


26.47 


68.91 .06619 


9 


39.33 .42495 


. 06561 


1.16650 


2.00 


38.76 .00500 


10 


69.29 .23339 


.03736 


1.21282 


30.15 


'60 .'29 .07539’ 


11 


58.63 .27522 


. 04414 


1 . 23549 


2.00 


57.70 .00500 


12 


48.89 .34594 


.05293 


1 . 22826 


2.00 


40.19 .00500 


13 


38.50 .41329 


. 06718 


1.19430 


2.0 0 


37.94 .00500 



FRONTAL AREA 


AFR 


.06264 


flow length 


L 


.16670 


MATRIX density 


(LS/FT3) 


70.51272 


• free Flow area 


AC 


.05040 



SOLID FRONTAL AREA AS .00424 



POROSITY “CR .84000 

HEaT~ TRANSFER AREA A'“ 10759300“ 

HYDR AU LIC R A D I US Rh .00046 

AREA DENSITY (A/VM) 0ETa' 1604*5 

COND.PARAMETE R ”C 0 ft RE CTl oN ' . 828 4 4“ 



*V- 
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TABLE VII. Summary of Heat Transfer and Friction Results 





.002 IN. SOLID NICKEL 


MOD. PARALLEL PLATE 






FRICTION DAiA 


HEAT TRANSFER DAIA 






CCNO 


MAX 




RUN ■ 


" NR " F' PARAM 


■ SLOPE NTO NR ■ ” 


J 



1 740.42 .02836 

2 608.95 .03389 

3 481.75 .04027 


.00925 
. 01125 
. 01420 


.57103 2.79 
.50706 3.00 
.61310 3.52 


729.70 

600.20 

474.04 


.00653 

.00722 

.00025 


4 373.92 .04731 

5 286.81 .05896 

6 187.14 .00669 


.01032 
. 02308 
. 03660 


' .’67679^ 4.63 

.71940 5.43 

.01039 7.74 


'360.55 

202.69 

104.45 


.01005 

.01271 

.01012 


7 115.43 .12779 

8 65.00 .21276 

9 3 ^.7 2 .3 72 3 3 


. 05931 
.10510 
.19163 


.06112 9.27 
.01110 7.70 
.77000 5.55 


113.77 

64.15 

35.21 


.02170 

.01023 

.01299 


FRONTAL AREA 


AFR 


. 06953 






flow length 


L 


. 16670 






MATRIX“DENSrTY 


aB/FT3) 


77.63067 






free flow area 


AC 


. 05939 






SOLID FRONTAL AREA 


AS 


.01014 






POROSITY ' ” 


PCR 


.05420 






HEAT transfer area 


A 


20.10750 






hydraulic radius 


RH 


.00049 
17A9 1 






AREA DENSITY (A/"VM) 


dE fA 








COND. PARAMETER CORRECTION 


1.00000 
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Maximum Slope 
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TABLE VIII. Ntu as a Function of Maximum Slope and Longitudinal 

Conduction Parameter 



APPENDIX A 



DATA REDUCTION RELATIONSHIPS 



GEOMETRY 



Porosity ■ Matrix free flow area 
Matrix frontal area 




Hydraulic diameter = 4x Free flow area 

Wetted perimeter 

Area density = Heat transfer surface area 

Matrix volume 





^ - Y»n " Apr L 

MASS RATE OF FLUID FLOW 



(A-1) 



(A-2) 



(A-3) 



The mass flow rate, m, is calculated in accordance with references 
[21] and [23]. 

m » 359 If c/o^ Fk Y V (A-4) 

Where: 



k = 

C = 

0 = 

d.= 

Y = 

Y = 
P = 
R = 
T = 



flow coefficient including velocity of approach 

orifice coefficient of discharge from [21] 

ratio of orifice diameter to inside pipe diameter 

orifice diameter 

thermal expansion factor 

expansion factor 

P 

the specific weight of the fluid flowing 

absolute static pressure at the orifice (Ibf/sq ft) 

53.35 ft IbF , the universal gas constant for air 
Ibm^R 

absolute temperature (®R) at the orifice 
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APo - Pressure drop across the orifice in inches of water 
Making substitutions for k and Y , (A-4) becomes; 

C do^Ra- Y 



m s 






r - PT 

In accordance with [23] figures (38) and (40b): 

Ya. » 1.0 

Y«i.o 

Also: 



(A-5) 



Where: 



and: 



P = ( Patti - if e '> ^ (Ibf/sq ft) 

is measured in inches of mercury 
Po is measured in inches of water 



T - t® + 459.7 

where is measured in degrees Fahrenheit 

Introducing the values and expressions above into (A-5) produces 
the working formula: 



m 



REYNOLD’S N13MBER 
By definition 



. C I p 



( Ao-t 45^.7) 



(A-6) 






Dh ^ 



Substituting for Dh * ^ 



u - 4Lm 

■ 7 m; 



(A-7) 
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FANNING FRICTION FACTOR 



— xim/7 //////////// ?\ — 

//)/f/nT777i 
n'r>}f ) / ? 7-7->-^fr>> t HA 



^ o 



The following equation describes the flow through 



th« natrix' [12J: 



R 



p. L ^ 






kc+Z-P’) 



ACCEUeR^TiChi 



f ^Vj£L 
Ac VT, 




where p = porosity ( 

It is assumed that and since the associated 

pressure changes are very small compared with the total pressure. It 

RT 

is also assumed that the perfect gas low applies ( ^ir= ). Making 

the necessary substitutions and solving for f produces for the iso- 
thermal case: 



f = [zPr, ^ - 




(A-8) 



Since the first term within the brackets above is much larger than 
the other terms, the following approximation can be made: 




Thus equation (A-8) becomes: 




2. Ptn '^2 



(kc + /<e) 




(A-9) 
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k c and /<a , the entrance and exit loss coefficients are functions of 
porosity, matrix configuration, and Re3mold*s number, and are obtained 
in accordance with reference [11] • 



MAXIMUM SLOPE 



The maximum slope of the generalized cooling curve as previously 
mentioned, is a function of Ntu and conduction parameter, X : 



Vtr, -t; / 

<1 C'W) 



“‘M4X 



(p(S)Luj\) 



(A-10) 



where: 



"c = -M— e 



WsC 






Therefore: 



and 

hn Cp 



and: 



'ty, „ mcp a 



d = ^de 

Cs 






fA-ll) 



also: 



and: 



fcPi 



tFt ~ti 

Fi “t; / 



•^1 



combining (A-11) and (A-12) yields: 

[tfj-t! 



“ 1 . J(^Fx-t:) 






- 









de 



MAX 









(A“ia> 



(A-13) 
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Now consider the following sketch of a generalized cooling curve: 




60 



MAX 




~ CHAftT SPeED 

d (tpj-tf.) = V 



Combining the above produces the maximum slope: 







-^MAK 



Os I V 

— K -^t^KCH^lLT ^P£EO 

Cr ^ ^ 



(A-14) 



Ntu can now be determined from conduction parameter X maximum 
slope. 
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COLBURN j 



Colburn's heat transfer parameter, j, is defined: 

j = Mst 

&Cp 

thus: 



(A-15) 



J = Htu ^ Npr*'* 



(A-16) 



From equations (A-1) and (A-3): 4 c 

pc I ■ Sha 

Apr 



)S= i- 

Af,/. 



J= Nt« -f- I 



(A-17) 



From equation (A-17) it can be seen: 

HEAT TRANSFER POWER 

The heat transfer power per unit area per degree temperature 
difference is [12]: 

Cp/* J_ 

Oh 



h = -TT^. -ir j Wr 



(A-18) 



For dry air at standard conditions equation (A- 18) becomes: 



llsro = .O^I95• M 



Nr o 



Dh 
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FLOW FRICTION POWER 



The flow friction power per unit area is defined as fl2]: 



E = 




(A-20) 



For dry air at standard conditions equation (A-20) becomes: 

,3 .3 



•STD 






s%.rt 



(A-21) 



For comparison purposes the surface geometries were reduced to a 
common hydraulic diameter. 
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APPENDIX B 



DESCRIPTION OF EQUIPMENT 

PRIME MOVER ; 

Air as the working fluid is drawn through the system by means of 
a 30 HP, multistage, Spencer Turbo-Compressor. This compressor is de- 
signed for a flow rate of 550 cubic feet per minute. 

FLOW METER ; 

Flow rates are measured by means of an ASME standard orifice section 
designed for convenient changing of orifice plates. Thin plate concen- 
tric orifices with throat diameters ranging from .308 to 2.31 inches are 
employed allowing for a wide range of flow rates. Pressure taps are 
located d upstream and d/2 downstream. Inside pipe diameter, d, is 
3.08 inches. 

HEATER SYSTEM ; 

The heater section is made up of 28 grid type nichrome wire heaters. 
Two heaters are wound to each bakellte frame and thus may be energized 
in pairs. A schematic wiring diagram for one heater frame is shown below: 






VOLTA6E 
500 RCE 



\f 





HeATcM. SVSTcM 
Cti - OFF 
SWiTCff 



( 



^CLccroiz 

SWITCH 



1 



I 



f 

f 







I 

1 

I 

i. 
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When energized, all heaters operate in parallel, the number of 
heaters in use depending upon flow rate. 

For maximum slope operation the voltage source is 220 volts AC and 
can be varied by means of a General Radio Company *Variac** Autotrans- 
former. 

For cyclic operation the voltage source is a 1.5 KW General Electric 
Amplidyne connected in series with a D.C. source. The resulting input 
to the heaters is a sine wave voltage with an amplitude of 50 volts and 
an average value of 50 volts. Frequencies of .1 to 1.0 cycles per second 
are employed with a much wider range available. 

MATRIX HOLDER MD TEST SECTION ; 

This section is constructed of polyethylene machined to close toler- 
ances. The T 2 and T^ thermocouple grids and the matrix upstream and 
downstream pressure taps are mounted in the frame. The T^ thermocouple 
grid is mounted in the removable matrix holder. sheet styrofoam in- 
serts provide insulation and a snug fit for the test matrix. The flow 

✓ 

channel measures 3-3/16** X 3-3/16** and matrices up to 3** in length may 
be accommodated in the test section. (See Figure 12) 

TEMPERATURE MEASUREMENT ; 

The temperature at the orifice meter is measured by means of a 
copper-constantan thermocouple and read in millivolts with a Rubicon 
Portable Precision Potentiometer. 

All other temperatures are measured by means of iron-constantan 
thermocouples, each temperature being measured by a group of five 
thermocouples connected in series. The T^ thermocouples are bound to- 
gether and insulated from each other by teflon tape. This bundle is 
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mounted at the throat of the intake cone, upstream of the wire screen 
flow straightener. The T^f T^jand grids are connected to a multiple 
barrel switch so that the emf generated by may oppose the emf generat- 
ed by either or Thus the difference between ambient and matrix 

inlet or exit temperature can be read* These temperature differentials 
are recorded on a Minneapolis-Honeywell "Brown** Strip Chart Recorder with 
variable sensitivity and chart speed. For maximum slope operation, 
versus T^ is recorded. For cyclic temperature variation, T^^ versus T 2 , 
and versus T^ are recorded. 

PRESSURE MEASUREMENT : 

The four pressure taps shown in Figure 8 are connected to Ellison 
draft gages by means of Imperial Company **Poly-flo** plastic tubing. The 
draft gages are a combination of inclined and vertical gages. They are 
shown in Figure 9. For small pressure differentials, Vernon Hill Company 
Type **C** Micromanometers are used. Pressures of 1.25** to 3.0** of water 
are read by means of the inclined draft gages. Pressures below this range 
are measured by means of the micromanometers; Pressure above 3.0** are 
read on vertical draft gages. 
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APPENDIX C 



CONDUCTION PARAMETER FOR PERFORATED MATERIAL 
The porous material used in the perforated plate matrices investi- 
gated in this report was made by Perforated Products, Inc. by means of 
an electro-depositing process. 

For a solid plate or fin, the conduction parameter is defined as: 

\ - 

^ LC 

where A is the solid cross-sectional area for conduction which is 
s 

identical to the product of plate width and thickness, and L is the 
conduction path length which is identical to the matrix flow length. 

For a perforated material such as the type 160/40 TV shown in 
Figure 38(a), the cross-sectional area for conduction varies in the 
flow direction, and the conduction path length is greater than the flow 
length . 

The definition of conduction parameter is based upon a solid materi- 



al, i . e . , 



V - _ _L_ 

LC " RthC 



where R . = — r the matrix thermal conduction resistance, 

th k A 
s s 

Using the above definition for a perforated material, an "equivalent** 
conduction parameter, X|^ , can be evaluated by determining the thermal 
resistance defining a relative shape factor S*, the conduc- 

tion parameter of a perforated material, X(c , can be related to the con- 
duction parameter of the same material as if it were unperforated: 

\ V C At, 

: X k = ^ A ^ 



or; 
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By the analogy between conduction heat flow and electric current, note 
that for similar shapes, the ratio of electrical resistances: 

Bjk « s' 

0^n)i< 

In Figure 38(b) are shown two shapes, each representing the shaded 
portion of Figure 38(a)* The shape on the left represents an unperfor- 
ated material, and the shape on the right is an incremental model for the 
160/40 TV material. The relative shape factor is determined from the 
ratio of the electrical resistances of these two models. Since the 
material is of uniform thickness, two dimensional models can be construct- 
ed from electrically conducting paper (such as Teledeltos paper) as shown 
in Figure 38tc) , with low resistance silver paint applied as indicated. 

The electrical resistances are then measured with an ohm-meter. This 
analog method produced the following results for the type 160/40 TV per- 



forated nickel material: 

^ s ,3/5- 

^ tu. 2.7 



8.5" 



thus: V 

Whereas, hk = .127X as determined in reference [22] using: 

Ak L 






where: 



Z Ak.1 

Z i-ki 

and: and are approximations indicated in Figures 38(a), and 

(b). 

The effects of these approximations upon Ntu and j are shown in 
Figure 15. 
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APPENDIX D 



EXPERIMENTAL UNCERTAINTY INTERVALS 

According to reference [13], the uncertainty interval for a 
result is: 



Wn '^2’’**'^n uncertainty intervals of the factors making 

up the result. 

R is the expression for the result 

V,, are the factors making up the result. 

1 z n 

The expression for Colburn j is: 




(D-l) 






The uncertainty interval for Colburn j is: 





For Ntu = 3.0; X= 0 etc: 
= 7.5% 




0 



For Ntu - 25.0; X ■ .05 etc: 




10 . 2 % 
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APPENDIX E 



DIGITAL COMPUTER PROGRAM FOR DATA REDUCTION 
Included herein is a computer program for a CDC 1604 digital com- 
puter. By means of this program all necessary calculations were made 
to produce the heat transfer and flow friction data appearing in the 
tables and graphs of this report. 

A data reduction program was written by Bannon [1]. This program 
was then further sophisticated by Piersall [22], producing conduction 
parameter X , maximum slope, and friction factor as well as other flow 
information from geometric and experimental data inputs. 

The present program will now produce all necessary heat transfer and 
flow friction data in table form from a minimxim of input. When data from 
the single-blow technique is processed, a cuirve-f itting interpolation 
subroutine is utilized to determine Ntu from maximum slope and conduction 
parameter. The input for this subroutine is Howard*s [9] maximum slope 
data (Table VIII). For the cyclic technique, Ntu is determined from 
equation (2-33) by means of an iteration process. 

The program will process data obtained by means of either the 
single-blow or cyclic technique and data from more than one matrix can be 
reduced in each computer run. 
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IF (FR-. 001)61, 61, 69 

69 EM=(D(M)*' .24)/(WM*CM*6.2832*FR*3600.0) 
CALL CYCLE (RB,EM,p) 

TU(M)=P 
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63 P = P + .002 
60 TO 62 

64 RETURN 
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